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Abstract 

O 

We examine the flavor dependence of the nuclear modification factor Raa i n the pQCD 
calculations at LHC energies. The computations are performed accounting for radiative 
and collisional parton energy loss with running coupling constant. Our results show 
that the recent LHC data on the Raa for charged hadrons, .D-mesons and non-photonic 
electrons agree reasonably with the pQCD picture of the parton energy loss with the 
rv 1 dominating contribution from the radiative mechanism. 
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1. The parton energy loss in the quark-gluon plasma (QGP) is widely believed to be a 
source for strong suppression of high-pr hadrons in AA-collisions (usually called the jet 
quenching) observed at RHIC and LHC. Understanding the underlying mechanisms of the 
parton energy loss is of great importance for application of the jet quenching to probing the 
hot QCD matter produced in AA-collisions. In the pQCD picture fast partons lose energy 
mostly due to induced gluon radiation [1, 2, 3, 4, 5, 6]. The effect of collisional energy loss 
[7] for the RHIC and LHC conditions is likely to be relatively small [8, 9]. Unfortunately, 
uncertainties in the pQCD-based models of the jet quenching remain large (mostly due to 
difficulties in modeling multiple gluon emission). For the nuclear modification factor Raa 
they are perhaps about a factor two. Despite this, it seems relatively safe to assume that 
predictions for variation of the Raa should be more robust, if the parameters are already 
adjusted to fit some set of experimental data. 

From the point of view of the underlying physics of the jet quenching it is very in- 
teresting to compare Raa for light and heavy flavors. It was suggested [10] that for the 
heavy quarks the dead cone effect should suppress induced gluon emission and give rise 
to an increase of the Raa- However, the observed at RHIC strong suppression of the non- 
photonic electrons from the B/D— meson decays [11, 12, 13] seemed to be in contradiction 
with this picture. It may indicate that for RHIC conditions the dead cone suppression 
is not very strong or that the radiative mechanism is not the dominating one at all. It 
stimulated the renewed interest in the collisional energy loss [14]. Although, by adjusting 
the coupling constant one can obtain a sufficiently strong heavy quark suppression due 
to the collisional mechanism alone, this scenario does not seem to be realistic (at least 
for pt > 5 — 10 GeV). Calculations of the radiative and collisional energy losses with the 
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same a s and the Debye screening mass performed in [8] clearly demonstrate that the col- 
lisional loss is relatively small for relativistic partons and unlikely to change significantly 
the heavy quark energy loss (see also [9]). 

In [10] the dead cone suppression was estimated from a qualitative analysis neglecting 
the quantum finite-size effects. Calculations of the induced gluon emission from heavy 
quarks in a brick of QGP [15] within the light-cone path integral (LCPI) approach [2], 
which treats accurately the mass effects, demonstrate that at energy ~ 10 — 20 GeV for 
c-quark the induced gluon spectrum is very similar to that for light quarks and AE C ps 
AE Uj d, s , and only for for 6-quark the gluon emission is suppressed (but not so strongly as 
predicted by the dead cone model [10]). At high energies (> 100 — 200 GeV) the radiative 
energy loss has an anomalous mass dependence with AE b > AE C > AE UjdjS due to the 
quantum finite-size effects in radiation of hard gluons [15]. In light of these results we can 
expect that the nuclear modification factor for the heavy quark jets for RHIC and LHC 
conditions should be qualitatively similar to that for light partons already at ~ 10 — 20 
GeV. Although accurate simulations and comparison with experiment are needed to reach 
definite conclusions. 

In the present work we examine the flavor dependence of the nuclear modification 
factor within the LCPI approach [2] and compare our results with the latest LHC data 
on the Raa for charged hadrons [16, 17], D-mesons [18, 19] and non-photonic electrons 
[20] in Pb + Pb collisions at yfs = 2.76 TeV. In evaluating the nuclear modification factor, 
besides the radiative energy loss, we include the collisional one. Both the radiative and 
collisional contributions are calculated with running a s . We account for accurately the 
fluctuations of the parton path lengths in the QGP. We find that the predicted flavor 
dependence of the Raa agrees reasonably with the LHC data. 

2. We calculate the nuclear modification factor employing the method developed in Ref. 
[21], to which the interested reader is referred for details. Here we just outline the main 
aspects of the calculations necessary for understanding of our strategy and interpretation 
of the results. 

For a given impact parameter b the Raa can be written as 

dN(A + A-+h + X)/dp T dy 
A ) T AA (b)da(N + N h + X) / dp T dy ' [ ) 

Here is the particle transverse momentum, y is rapidity (we consider the central region 
y — 0), TaaQ)) = J dpT A (p)T A (p — b), T A is the nucleus profile function. The differential 
yield in AA-collision can be written in the form 

dN(A + A^h + X)_ l ipTA(p)TA(p _ h) d* m (N + N->h + X) 



dp T dy J dp T dy 

da m (N + N -^h + X) f x dz . n ,da(N + N^i + X) 



dprdy ^ Jo z 2 dp l T dy 

Here p l T = Pt/z is the parton transverse momentum, da(N + iV — >■ % + X)/dp l T dy is the 
hard cross section, D™^ is the medium-modified fragmentation function (FF) for transition 
of a parton % into the observed particle h. For the parton virtuality scale Q we take the 
parton transverse momentum p l T . 
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We assume that the induced radiation stage occurs after the DGLAP stage which gives 
the input parton distribution for the induced gluon emission stage. It seems reasonable 
since for jets with E < 100 GeV the typical time scale for the DGLAP stage is relatively 
small (< 0.3 — 1 fm [21]), and in first approximation it is legitimate to neglect interference 
of the DGLAP and the induced gluon emission stages. Symbolically the medium-modified 
FF reads 

D h/iiQ) « D h/j (Q ) ® Df /k ® D k/l {Q) , (4) 

where ® denotes ^-convolution, D k /i is the ordinary DGLAP FF for i — > k parton tran- 
sition, DJJ k is the FF for j — >■ k parton transition in the QGP due to induced gluon 
emission, and -D/i/j describes parton hadronization outside of the QGP 1 . In (4) Q is the 
scale at which the DGLAP parton showering is stopped. As in [21] we take Q = 2 GeV. 

We computed the DGLAP FFs with the help of the PYTHIA event generator [24]. 
The one gluon induced spectrum, dP/dx, was calculated within the LCPI approach [2] 
employing the method developed in [25]. The DJ) k has been obtained from dP/dx ac- 
counting for multiple gluon emission within Landau's method as in [26]. Note that we 
include the q — >■ g transition as well, which is usually neglected. For the Dh/j(Qo) we use 
the KKP [27] FFs for light partons, and Peterson FF for heavy quarks (with parameters 
e c = 0.06 and e b = 0.006). For the non-photonic electrons we evaluated the FFs c — > e 
and b — > e treating them as the two-step fragmentations c — > D — > e and b — > B — > e. 
The distributions B/D — >■ e were calculated using the CLEO data [28, 29] on the electron 
spectra in the S/D-meson decays. We neglected the B — >■ D — >■ e process, which gives a 
negligible contribution [30]. 

The hard cross sections were calculated using the LO pQCD formula with the CTEQ6 
[31] parton distribution functions. To simulate the higher order effects we take for the 
virtuality scale in a s the value cQ with c = 0.265 as in the PYTHIA event generator 
[24]. This prescription allows us to reproduce well the pr-dependence of the spectra in 
pp-collisions 2 . In calculating the Raa we account for the nuclear modification of the 
parton densities (which leads to some small deviation of Raa from unity even without 
parton energy loss) with the help of the EKS98 correction [32]. 

As in [21] we take m q = 300 and m g = 400 MeV for the light quark and gluon 
quasiparticle masses supported by the analysis of the lattice data [33]. For the heavy 
quarks we take m c = 1.2 GeV and nib = 4.75 GeV. We use the Debye mass obtained in 
the lattice calculations [34] that give the ratio //d/T slowly decreasing with T (hd/T w 3 
at T ~ 1.5T C , fi D /T w 2.4 at T ~ 4T C ). 

We use the running a s frozen at some value a{ r at low momenta (the technical details 
for incorporating the running a s can be found in [25]). For gluon emission in vacuum a 
reasonable choice is a{ r ~ 0.7 [35, 36]. The RHIC data on the pion Raa in Au + Au 
collisions at = 200 GeV support a{ r ~ 0.5 — 0.6 [21]. But the analysis [37] of the first 

1 Thc approximation (4) ignores creation in the QGP of the anomalous jet color states, which may be 
important for the baryon Raa [22] at not very high px and the jet structure in the soft region [23]. But 
it should be reasonable for evaluating the Raa for charged hadrons, which is dominated by the charged 
pions, and the Raa for heavy flavors. 

2 Although we use the LO formula for the heavy quark cross sections, the p-r-dependences (and the 
c/b ratio) of our cross sections agree well with the more sophisticated FONLL calculations [30] (the 
normalization of the cross sections is unimportant for the Raa at all). 
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LHC data on the Raa for charged hadrons in Pb+Pb collisions at y/s = 2.76 TeV obtained 
by ALICE [38] shows that they agree better with a{ r ~ 0.4 — 0.5. The calculations with 
a fixed coupling constant [39, 40] also indicate that it can be smaller at LHC energies. 
The reduction of a{ r from RHIC to LHC is probably a manifestation of the thermal 
suppression of a s due to the growth of the initial temperature of the QGP at LHC. We 
will see that the new data from CMS [16] and ALICE [17] also support a{ r ~ 0.4 - 0.5. 

We view the collisional energy loss as a perturbation [21], and account for its effect 
simply by redefining the initial QGP temperature in calculating the radiative FF according 
to the condition 

AE rad (To) = AE rad (T ) + AE col (T ) , (5) 

where AE rad / co i is the radiative/collisional energy loss, T is the real initial temperature of 
the QGP, and T is the renormalized temperature. We solve (5) in linear approximation in 
T ' 3 -T 3 , which gives T ' 3 = T 3 + AE col (T )/[dE rad (T )/dT$]. It was done for each parton 
trajectory in the QGP (separately for quarks and gluons). The collisional energy loss has 
been evaluated in the Bjorken method [7] with an accurate treatment of kinematics of the 
binary collisions (the details can be found in [8]). 

3. We perform the computations for Bjorken 1+1D longitudinal expansion of the QGP 
[41], which gives TqT = T 3 r. We take r = 0.5 fm. For simplicity we neglect variation 
of the initial temperature T in the transverse directions. We evaluated T using the 
data on the charged hadron multiplicity pseudorapidity density dN ch /dr] [42, 43] and the 
entropy/multiplicity ratio dS / dy j dN c h/ dr] 7.67 obtained in [44]. It gives To ~ 420 
MeV for central Pb + Pb collisions at = 2.76 TeV. For each jet we calculate accurately 
the fast parton path length in the QGP, L. To take into account the fact that at times 
about 1—2 units of the nucleus radius the QGP should cool quickly due to transverse 
expansion [41], we impose the condition L < L max . We performed the computations for 
Lmax = 8. The bigger value L max = 10 fm gives almost the same. 

In Fig. 1 we compare the theoretical Raa for charged hadrons obtained for a{, r = 0.5 
and 0.4 to the data from CMS [16] and ALICE [17] for 0-5% central Pb + Pb collisions at 
y/s = 2.76 TeV. The results are presented for the radiative mechanism alone and with the 
collisional energy loss. We show our results for p T > 5 GeV since at smaller momenta our 
perturbative treatment is hardly applicable. Fig. 1 shows that the collisional mechanism 
suppresses Raa by ~ 20% at pt ~ 10 GeV, and ~ 10% at pt ~ 100 GeV. One sees that the 
teoretical Raa (for radiative plus collisional energy loss) for the window a{ r ~ 0.4 — 0.5 
agrees reasonably with the experimental data. The agreement is somewhat better for 
a{ r = 0.4. 

In Fig. 2 we compare our results with the ALICE data [18, 19] on the Raa for D- 
mesons in Pb + Pb collisions at \fs = 2.76 TeV for 0-20% and 0-7.5% centrality bins. 
Fig. 2 shows the results for the c — > D fragmentation. We found that the effect of the 
6-quark (due to b — >■ B — )■ D vacuum fragmentation) increases the Raa only by about 
2%. From Fig. 2 we can conlude that the same window in a{ r as for light hadrons allows 
to obtaind a fairly reasonable description of the D-meson data as well. 

In Fig. 3 we compare our calculations of the Raa for non-photonic electrons with 
the recent ALICE measurement [20]. We show the contibution from the charm and 
bottom quarks separately and the total electron Raa- Note that for the bottom quark 
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our treatment of the collisional mechanism as a pertubation to the radiative one, with 
the help of (5), loses accuracy at < 5 — 6 GeV. In this region the collisional correction 
becomes too large for the predictions to be robust. It happens since the Raa becomes 
sensitive to the low energy region where for the bottom quark AE co i > AE ra d. Evidently, 
in this regime the radiative and collisional mechanisms must be treated on an even footing. 
Unfortunately, this problem remains unsolved. For the charm quark this complication does 
not arise since accross the whole energy range the collisional energy loss remains relatively 
small [8]. Fig. 3 shows that at ^,6 — 7 GeV our results agree with the data fairly 
well. Note that for the RHIC conditions our results also agree reasonably with the data. 
For the 0-5% central Au + Au collisions at yfs = 200 GeV for a f / ~ 0.5 - 0.6 (what is 
needed for agreement with the Raa for pions) at pr ~ 6 — 8 GeV our calculations give the 
electron Raa ~ 0.25 — 0.35, which agrees reasonably with the STAR [12] measurement. 
A detailed discusion of the non-photonic electrons for RHIC and LHC energies will be 
given in a forthcoming publication. 

Thus, our pQCD model with the radiative energy loss combined with relatively small 
collisional energy loss gives a reasonable description of the latest LHC data on the Raa 
both for light and heavy flavors at px > 5 GeV. 

4. In summary, we have examined the flavor dependence of the nuclear modification 
factor Raa in the pQCD picture and checked its consistency with that observed at LHC. 
We show that the LHC data on the Raa for charged hadrons [16, 17] and D-mesons 
[18, 19] in central Pb + Pb collisions at y/s = 2.76 TeV can be reasonably described in 
the pQCD scheme, universal for light and heavy flavors with relatively small collisional 
energy loss. We found that the ALICE data [20] on the Raa for non-photonic electrons 
can be described fairly well in our model as well. 

We conclude that the recent LHC data on the Raa for the light and heavy jets give 
strong support for the validity of the pQCD parton mass dependence of the energy loss 
with relatively small effect of the collisional mechanism. The collisional mechanism be- 
comes very important only for the bottom quark at momenta < 6 — 8 GeV. For accurate 
pQCD calculations in this region a better understanding of the interplay of the radiative 
and collisional mechanism is required. 
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Figure 1: The nuclear modification factor for charged hadrons at y = for 0-5% central 
Pb + Pb collisions at y/s = 2.76 TeV for a{ r = 0.4 (upper curves) and 0.5 (lower curves). 
The solid line shows the calculations with the radiative and collisional energy loss, and 
the dashed line shows the results for the radiative mechanism alone. The experimental 
points are the data from CMS [16] (circles) and ALICE [17] (squares). Systematic errors 
are shown as shaded areas. 
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Figure 2: The D-meson nuclear modification factor for 0-20% (left) and 0-7.5% (right) 
central Pb + Pb collisions at y/s = 2.76 TeV for a{ r = 0.4 (upper curves) and 0.5 (lower 
curves) at y — 0. The solid line shows the calculations with the radiative and collisional 
energy loss, and the dashed line shows the results for the radiative mechanism alone. The 
experimental points are the ALICE data [18] (left panel), [19] (right panel) for average 
D°,D + ,D* + . Systematic errors are shown as shaded areas. 
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Pb+Pb, s 1/2 =2.76 TeV, 
+ ALICE preliminary 
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Figure 3: The electron nuclear modification factor for 0-10% central Pb + Pb collisions 
at -y/s = 2.76 TeV for a{ r = 0.4 (upper curves) and 0.5 (lower curves) at y — 0. The 
solid line shows the total Raa, the dotted and dashed lines show the Raa for charm and 
bottom contributions, respectively. The experimental points are the preliminary ALICE 
data [20]. Systematic errors are shown as shaded areas. 
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